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Abstract 

In this paper we introduce a new flexible class of distributions with bounded support, called re¬ 
flected Generalized Topp-Leone Power Series (rGTL-PS), obtained by compounding the reflected 
Generalized Topp-Leone (van Drop and Kotz, 2006) and the family of Power Series distributions. 
The proposed class includes, as special cases, some new distributions with limited support such as 
the rGTL-Logarithmic, the rGTL-Geometric, the rGTL-Poisson and rGTL-Binomial. This work 
is an attempt to partially fill a gap regarding the presence, in the literature, of continuous dis¬ 
tributions with bounded support, which instead appear to be very useful in many real contexts, 
included the reliability. Some properties of the class, including moments, hazard rate and quantile 
are investigated. Moreover, the maximum likelihood estimators of the parameters are examined 
and the observed Fisher information matrix provided. Finally, in order to show the usefulness of 
the new class, some applications to real data are reported. 
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1 Introduction 

In recent years, many authors have focused their attention on the proposition of new and more flexible 
distribution functions, constructed using various transformation techniques such as, for example, the 
Beta-generated distribution by Eugene et al. (2002) and Jones (2004), the Gamma-generated distri¬ 
bution by Zografos and Balakrishnan (2009), Kumaraswamy-generated distribution by Cordeiro and 
de Castro (2011), McDonald-generated distribution by Alexander et al. (2012), Ristic and Balakrish¬ 
nan (2012), Weibull-generated distribution by Bourguignon et al. (2014), just to name a few. Such 
transformation techniques may be viewed as special cases of the Transformed-Transformer method 
proposed by Alzaatreh et al. (2013). Other proposals are based on the Azzalini’s method (Azzalini, 
1985, 1986) and its extensions (Domma et al., 2015). Using these techniques, many of the known 
distributions such as, for example. Normal, Exponential, Weibull, Logistics, Pareto, Dagum, Singh- 
Maddala, etc., have been generalized. 

Almost all these new proposals concern distributions with unbounded support. In the face of the nu¬ 
merous proposals of distributions with unbounded support emerges, undoubtedly, the great scarcity 
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of distributions with bounded support (Marshall and OUcin, 2007, pag. 473), although there are many 
real-life situations in which the observations clearly can take values only in a limited range, such as 
percentages, proportions or fractions. Papke and Wooldridge (1996) claim that variables bounded 
between zero and one arise naturally in many economic setting; e.g. the fraction of total weekly hours 
spent working, the proportion of income spent on non-durable consumption, pension plan participa¬ 
tion rates, industry market shares, television rating, fraction of land area allocate to agriculture, etc. 
Various examples of proportions in the unit interval used in empirical finance are discussed in Cook 
et al. (2008). 

Also in reliability analysis, different authors refer to continuous models with finite support in order 
to describe lifetime data. This is often motivated by considering physical reasons such as the finite 
lifetime of a component or the bounded signals occurring in industrial systems (see, for example, 
Jiang, 2013; Dedecius and Ettler, 2013). In this perspective, the models with infinite support can be 
viewed as an approximation of the realty. Furthermore, when the realiability is measured as percent¬ 
age or ratio, it is important to have models defined on the unit interval (Gen 9 , 2013) in order to have 
plausible results. 

It is well known that the most used distribution to model continuous variables in the unit interval 
is the Beta distribution. The popularity of this distribution is certainly due to the great flexibility of its 
density function, in fact it can take different forms such as constant, increasing, decreasing, unimodal 
and uniantimodal depending on the values of its parameters. This distribution has been used in var¬ 
ious fields of science such as, for example, biology, ecology, engineering, economics, demography, 
finance, etc. (for a detailed discussion see Johnson et al. (1995), Nadarajah and Kotz (2007)). On the 
other hand, mathematical difficulties underlying the use of this model, due to the fact that its distribu¬ 
tion function cannot be expressed in closed form and its determination involves the incomplete beta 
function ratio, are well known. 

Recently, several authors have proposed an alternative to the Beta distribution by recovering the dis¬ 
tribution proposed by Kumaraswamy in 1980 in the context of hydrology studies. As pointed out by 
Jones (2009), Kumaraswamy’s distribution has many of the properties of the Beta distribution and 
some advantages in terms of tractability, in particular its distribution function has a closed form and 
it does not involve any special function. In fact, this distribution turns out to be a special case of the 
generalized Beta of the first type proposed by McDonald (1984) (see Nadarajah (2008)). 

Perhaps thanks to a work of Nadarajah and Kotz (2003), a renewed interest has been recently 
developed for another distribution defined on bounded support, the Topp-Leone (TL) distribution. 
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proposed by Topp and Leone (1955) and successively studied by different authors (see, for example, 
Ghitany, 2007; Gen 9 , 2012; Vicari et al., 2008) also in the reliability context (Ghitany et ah, 2005; 
Gen 9 , 2013; Condino et ah, 2014). 

Starting from a generalized version of the TL distribution, van Drop and Kotz (2006) proposed 
the reflected Generalized Topp-Leone (rGTL) distribution. Similarly to the Beta distribution, this 
distribution has a density that can be constant, increasing, decreasing, unimodal and uniantimodal, 
depending on the values of its parameters. Moreover, it has a strictly positive density value at its 
lower bound and a closed form of its distribution function. However, the hazard function of the rGTL 
shows a certain rigidity since it is always increasing. This fact is a real weakness of the model, in 
particular in the field of reliability theory and survival analysis. 

With the aim to propose a new flexible model defined on the unit interval, in this paper we consider 
the standard rGTL distribution and introduce the reflected Generalized Topp-Leone Power Series 
(rGTL-PS) class of distributions, obtained by compounding the Power Series distributions and the 
standard rGTL. 

In the recent literature, many new distributions are obtained by compounding a continuous distri¬ 
bution with a discrete one. An interesting motivation of this procedure can be found in the process 
underlying the failure of a series system composed by Z component. If F* is the lifetime for the i — th 
component, the system will fail when the first component fails, so the lifetime of the whole system 
is F = min(Fi,..., Yz). In this situation, assuming that the lifetimes of the components are indepen¬ 
dent, it is easy to obtain the probability of failure for the system by compounding the distribution of 
F’s and the distribution of Z, as follows: 

P{Y<y) = P[min(Fi,...,Fz)<|/] = 

OO 

= 1 - 5^P(Fi > y)P{Y 2 > y)...P{Yr, > y)P{Z = z). 

z=Q 

By considering different cumulative distribution functions for Z and for Yi, various models are pro¬ 
posed during the last years. It is the case, for example, of the exponential-geometric distribution 
(Adamidis and Loukas, 1998), of the exponential-Poisson-Lindley distribution (Barreto-Souza and 
Bakouch, 2013) and of the exponential-logarithmic distribution (Tahmasbi and Rezaeib, 2008), to 
name a few. Other distributions, such as the Weibull Power Series (Morais and Barreto-Souza, 2011), 
are obtained by describing the random variable Z through the Power Series distribution, or by con¬ 
sidering a similar procedure to that just mentioned, involving the maximum rather than the minimum 
of the lifetimes (Nadarajah et al., 2013; Flores et al., 2013). 

The rGTL-PS distribution obtained preserves the main advantage of the rGTL distribution with 
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respect to the Beta distribution, that is its cumulative distribution function has a closed form and 
therefore the quantile functions are easily obtainable and one can easily generate a random variable 
from rGTL-PS distribution. Furthermore, the rGTL-PS density function has an analogous flexibil¬ 
ity of the Beta density function, i.e. the shape of density can be increasing, decreasing, unimodal 
and uniantimodal. Finally, the shape of the hazard, besides being increasing and bathtub as the Beta 
distribution, also shows a N-shape (very useful in the context of reliability theory, see for example, 
Bebbington et al. (2009), Lai and Izadi (2012)). These properties characterize our proposal as a valid 
alternative to the Beta distribution. 

This article is organized as follows. In Section 2, we define the rGTL-PS distribution. Some proper¬ 
ties, such as the moments, the hazard rate and the quantile are derived in Section 3. The maximum 
likelihood estimation is discussed in Section 4 and some special cases are studied in Section 5. Finally, 
various applications on real data sets are reported in Section 6. 

2 Reflected Generalized Topp-Leone Power Series distribution 

A random variable Y is said to have a standard reflected Generalized Topp-Leone (rGTL) distribu¬ 
tion (van Drop and Kotz (2006)) if its cumulative distribution function {cdf) and the corresponding 
probability density function ipdf) are respectively given by 

G{y] a, z/) = 1 - (1 - yf [a - {a - 1)(1 - y)^ (1) 

and 

9 {y,o>,iy) = z/(l - yY~^[a -{a - 1)(1 -y)]’'~^[a - 2{a - 1)(1 - ?/)] (2) 

with 0 < a < 2,u > 0 and 0 < y < 1. The density of the rGTL can be strictly decreasing, strictly 
increasing or may possess a mode or an anti-mode, according to the values of the parameters. 

In order to define the new distribution, we consider a sequence of independent and identically 
distributed continuous random variables Yi, Y 2 , ■■■Yz, with distribution function G{y), where Z is a 
discrete random variable following a Power Series distribution truncated at zero, with probability 
function (p/) given by 

m = = z = G2,..- e>0 (3) 

where A{6) = finite and > 0. In Table 1 the most common distributions belonging 

to the Power Series family are reported. 
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Distribution 

Qz 

A{e) 

A{e) 

range of 0 

Logarithmic 

l/z 

- log(l - 6) 

(1-0)-^ 

(0.1) 

Geometric 

1 

1 

1 

(l-0)-2 

(0.1) 

Poisson 

llz\ 

— 1 

e® 

(0,oo) 

Binomial 

(T) 

(0 + 1)™ - 1 

m(0 + 1)™-^ 

(O.oo) 


Table 1: Some special cases of the Power Series distribution. 


Let V = min{Yi, ...Vz}- The conditional pdf of Y given that Z = z, i.e. Y\Z = 2 :, can be 
obtained from the distribution of the minimum of 2 ; random variables: 

f(y\Z^z) = z-[l-G(y)Y-^-g(y). (4) 

Hence, the joint pdf of (F, Z) is given by 

f{y, z) = f{y\Z = z)P{Z = z) = j^z ■ [1 - ■ g(y). (5) 

By denoting with A'(-) the derivative of A(-) with respect to the argument, the marginal pdf of Y 
is 


f\y) = 

Z = 1 


wf 


■ [1 - G{y)r^ 


■aiy) 


9g{y)A'{e[l-G{y)]} 

m 


( 6 ) 


with cdf given by 




(7) 


^^A{e) ' A{e) 

Replacing the expressions (1) and (2) in (6) and (7), we obtain the rGTL-PS pdf and the corre 
spending cumulative cdf as follows: 




f{y, tt, yS) = Y - i)(i - y)] (i - yT" Ma - (« - i)(i - y)] 

2=1 ^ ' 


PZ—l 


^ n 

E CtzG . 


IZZ] 


( 8 ) 


<7 f)^ 

2=1 ' ' 

°° „ Qz °° „ nz 

= ^-Yj^i^-G{y;a,uz)]=Yj^G{y,a,uz). (9) 

It is evident, from (8), that the density of the rGTL-PS class can be expressed as mixture of rGTL 
densities, g{y] a, vz), with weights In the following, we refer to this property because it enables 
us to obtain some mathematical properties of the rGTL-PS distributions, such as the moments. 
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3 Statistical Properties 


In this section, we study some properties of the rGTL-PS distribution. In particular, we determine the 
rth incomplete moment and ordinary moments using the known properties of the mixture distribu¬ 
tions. Moreover, we compute the quantile and the hazard rate, evaluating its behaviour at the extrems 
of the support. 

3.1 Moments 

In order to calculate the rth moment of rGTL-PS distribution, first we determine the incomplete 
moment of order r for a rGTL distribution then, by using the properties of the mixture distributions, 
we calculate the rth moment of a rGTL-PS distribution. 

Lemma 1 If Y ~ rGTL{a, v) then the incomplete moment of order r is 


E{Y^\Y<y*} = + 


-ky 


^T{u-h)h\ 


h=0 


{(2 — oi)By* (r + 1, z/ + /i) + 2 (q: — l)By* (r + 2, z/ + /i)} 


where y* G (0,1] and By* {a, b) = y°- ^(1 — yY ^dy. 

textbf/’roo/ Using (1 — kzY~^ = ^ ^ write 


E {Y'^\Y <y*} = ^ (2 — a) / w^{l — wY ^ [1 — k{l — w)Y ^ dw 


-|-2(q; — 1) / w^~^^{l — wY ^ [I — k{l — w)Y ^ dw 

Jo 


= oa 




/f=0 


vY-h)h\ 


+2(a-l) / w''^^{l-wY^'"~^dw 

Jo 

Corollary 2 IfY ~ rGTL{a, v) then the moment of order r is 


E{y'} = a''-‘r(t. + i)y; 




k)” r,o , 2(a-l)(r+l) 
- “> + (.. + /> +r+l)J 


B{r1,0h). (10) 


Proof It is enough to put y* = 1 in the Lemma. 
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From (10), it is easy to obtain the moment of order r of the rGTL-PS distribution: 


= Il^frGTL{Y^-,a,Uz) = 


Z=1 ' ' 


/ X 2 a-l r + 1 

-^TTTT 2-a +-^-^^ B(r + l,uz + h). 

r{uz-h)h\ ^ ^ {uz + h + r + 1) ^ ^ 


3.2 Hazard rate 


In this section, we verify that the hazard rate of rGTL-PS distribution is more flexible than the hazard 
rate of the rGTL. 

First of all, it is easy to show that the hazard rate of rGTL distribution is always increasing. Indeed, 
by (1) and (2) we obtain the hazard rate for the rGTL, as follows: 

, , , a — 2(a — 1)(1 — y) 

KaTLivic,,.) = ( 11 ) 

and, after easy algebra, the derivative of ( 11 ) with respect to y is given by 

dhrGTLiy-,0:,iy) V . 12 , / 1X2/1 x2\ 

- ay - +(«-!) ( 1 -!/)} 

that is positive Vy, u > 0 and a G (0, 2]. 

Starting from ( 6 ) and (7), we obtain the expression for the hazard rate of the rGTL-PS distribution. 


{[(a - 1)(1 -y)- + (a - 1)^(1 - yf] 


as follows: 


h{y;a,u,9) = 


d ■ g{y] a, u)A'{9[l - G{y; a, u)]} 


A{9[l-G{y,a,u)]] ' ^ ^ 

The hazard rate is a complex function of y. However, from (12) we have that for ?/ —)■ 0, the hazard 
rate tends to z/6*(2 — o;) ■ A'[9)/A{9), while for ?/ —)■ 1 the hazard rate tends to +cx). Indeed, we have: 

A'{9\i — G{y]a,i')]] 

\imh{y,a,iz,9) = Imi —— -—- ^9g{y,a,iz) 

A\9\i — G[ii',a,v)\\ 

= ai0 lim ——-— (13) 

y^i A{9[1 -G{y,a,v)]] 

where A'{9) = has the same radius of convergence of ^(6^) and x4'(0) = oi > 0. Observed that 
yl(0) = 0, if z/ < 1 then limj^_^i h{y, a, u, 9) = +cx), given that limj^_^i g{y; a, u) = +cx); while, if 
u > 1, given that limy_^i g{y; a, u) = 0, we have limj^_^i h{y; a, u, 9) = Now, using FHopital’s 
Rule we obtain: 

lim h{y, a, u, 9) = ai9 lim - 7 - 7777 ;- .7 (?/i ^ 

J/^I y^i yl'{6»[l - G{y, a, iy)]}[-9g{y, a, u)] 

r 9 '{y\a,v) 

= - lim -r. (14) 

y^i g{y;a,iy) 
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It can be shown that = —oo, so that limy_ 5 .i h{y; a, u, 9) = +cx). 

As we can see from Fig. 1, the hazard rate for the models belonging to the rGTL-PS class of 
distributions is much more flexible than the hazard rate of the rGTL distribution. In particular, besides 
having the increasing shape, it is also possible to have the bathtub shape and the upside-down bathtub 
and then bathtub shape. This wide range of different behaviours of the hazard rate makes the models 
belonging to the rGTL-PS class suitable models for reliability theory and survival analysis in cases of 
bounded domain. 


3.3 Quantile 

An advantage of the rGTL-PS distribution is the possibility to get the expression for the quantile func¬ 
tion. Indeed, the cdf for the rGTL-PS can be expressed in a closed form, as it can be noted from the 
expression (7) and therefore the quantile can be easily obtained by remembering the expression for 
the quantile of rGTL distribution: 


1 «- V a2-4(a-l)(l-p)l/-' 1 ^ ^ O 

^ 2(q-i) i ^ a ^ z 

yp= a = l 

ct+\/ 4 ( 0 —1)(1—p)!/" 

2(T^) 

A{e[l-G(yq\a,u)]} 


0 < a < 1 


and by putting q = F{yq-, a, z/, 6») = 1 

Thus, the quantile yq of the rGTL-PS distribution is given by 

Uq I 1 ^ , Q:, z/ I 


■l-^— l<a<2 


= < 


2{a-l) 

e 


a = 1 


1 -^^— 0 < a < 1 . 


2{a-l) 


4 Inference 


(15) 


(16) 


In order to estimate the parameters r] = {a, 12 , 9 ) of the rGTL-PS distribution, we consider the maxi¬ 
mum likelihood (ML) method. Lety = ( 1 / 1 , 1 / 2 , •••, l/n) be a random sample of size n from the rGLL-LS 
given by (6). The log-likelihood function for the vector of parameters y = (a, z/, 6*) can be expressed 
as: 

n n 

YVi y) = n\n9 — n\n A{9) + E ct, z/)-f In A [1 — G(|/i; a, z/)]} . (17) 

i=l i=l 
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Differentiating i{ri; y) with respect to a, u and 6, respectively, and setting the results equal to zero, 
we have: 


d^iviy) 

da 

d^{v;y) 

du 

dt{v,y) 

de 


ga(yi-,a,v) _ n S^n A {S[l-G(yi;a,i^)]} ^ / 7 /] — f) 

Z^i=l g(yy,a,iy) ^ ^i=l A'{0ll-G(yi;a,7y)]} ^ ^ 

_ g^{yi-a,u) _ n sr^n A {0[1-G(yi;a,iy)]} ^ / jA — f) 

Z^i=l g(yi;a,i.) ^ ^i=l A'{dll-G(yi;a,7y)]} ^ ^ 

n _ y.n A''{e[l-Giyy,a,u)]} m . \l q 

B MB') ^ l^i=l A' iBU-G( 7 ,.-a 7 ,\n iy)\ u 


^-G„{yi;a,u) = 0 


A {B[l-G{yi\a,v)W 


where the quantities g{.) and G{.) are reported in the Appendix. The system does not admit any ex¬ 
plicit solution; therefore, the ML estimates rj = {a,u,9) can only be obtained by means of numerical 
procedures. Under the usual regularity conditions, the known asymptotic properties of the maxi¬ 
mum l ik elihood method ensure that y/nifin — f]) A- iV(0, Sr,), where Sr, = [lim,r_roo T(? 7 )/n]“^ 
is the asymptotic variance-covariance matrix and I ( 77 ) is the Fisher Information matrix. Moreover, 
the asymptotic variance-covariance matrix of f) can be approximated by the inverse of the observed 
information matrix J{ri), whose entries are given in the Appendix. In order to build the confidence 
intervals and hypothesis tests, we use the fact that the asymptotic distribution of i) can be approxi¬ 
mated by the multivariate normal distribution, [J(i 7 )]~^), where [1(7?)]”^ is the inverse of the 

observed information matrix evaluated in i). 


4.1 EM algorithm 

In the original formulation due to Dempester et al. (1977), EM algorithm is a method for computing 
maximum likelihood estimate iteratively, starting from some initial guess, when the data are incom¬ 
plete. The EM algorithm has become a popular tool in the statistical estimation problems involving 
incomplete data, or in problems which can be posed in similar form. Each iteration of the EM al¬ 
gorithm consists of an Expectation (E) step, in which we calculate conditional expectation of the 
complete-data log-likelihood function given observed data, and a Maximization (M) step, in which 
we maximize this equation. 

In our case, we suppose that the complete-data Wi = (F,, Z,) (i = 1, ...n) consist of an observable 
part Yi and an unobservable part Zi. By (5) the corresponding log-likelihood function 

n n n 

oc \n{e)^Zi - \nA{d) + '^{zi - 1 ) ln[l - G{yi]a,u)] + '^lng{yi-,a,u). (19) 

2=1 2=1 2=1 

Then, the E-step of the algorithm requires the computation of the conditional expectation 

E[r(r/(^+i);W)|Y = y,r/(^)] (20) 

where is the current estimate of rj in the rth iteration. It is straightforward to 

verify that the E-step of an EM cycle requires the computation of the conditional expectation of the 
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conditional random variable (ZjV = y, Using (5) and (6), the conditional probability function 
of {Z\Y = |/,r/W) is 

9iz\y;r])= , . -^177-(21) 

A - G(y;a,iy)l} 

forzeR and y G (0,1). After easy algebra, the conditional expectation is given by 

E[Z\Y = V'->1 ^ _ Gfa + 1- (22) 

The M-step of EM algorithm requires the maximization of the complete-data likelihood over rj, with 
the unobservable part Zi replaced by their conditional expectations. Thus the at (r+l)th of EM 
is the numerical solution of the following nonlinear system: 


ar(??W;W) 

da 

9r(77W;W) 

du 

9r(77W;W) 

de 


2^1=1\Y ) ^1=1 \^l-g^[y-a('^) 

_ V” ( IN Gv{yi\a^'''> 9L^{yr,A^\’z(0'j _ 

0(r) l^i=l A{BM) ~ 


5 Special cases 

In this section, we furnish some results about special cases of the rGTL-PS class of distributions. In 
particular, by considering the quantities reported in Table 1 and the expressions given in (6) and (7), 
we obtain the pdf and the cdf for rGrL-Logarithmic, rGTL-Geometric, the rGTL-Poisson and the 
rGrL-Binomial distributions. Some plots of density functions and hazard rates are given in Fig. 1 to 
show the flexibility of these models. 


5.1 rGrL-Logarithmic distribution 

The pdf and the cdf of the rGrL-Logarithmic (rGTL-Log) random variable, respectively, are: 


f{y;a,u,e) = 


eujl - yy-Ya -{a- 1)(1 - y^-^c^ - 2{a - 1)(1 - y)] 
- log(l - e){l - e{l - yY[a - (a - 1)(1 - y)Y} 


F{y;a,u,9) = 1 


log{l - e{l - yY[a -(a- 1)(1 - y)Y} 


' log(l-0) 

with 0 < a <2,iy > 0 and 0 < 6^ < 1. The hazard rate is given by 


h{y;a,u,9) = 


_ dg{y,a,v) _ 

[1 - 9G{y] a, z/)] ■ [log(l - 9G{y] a, v)) ' 
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Some plots of the;? J/'and hazard rate are given in Fig. 1. Finally, considering that A ^[(1 —(;)74(6*)] = 
1 — (1 — 6Y~^, from (16), we obtain the quantile for rGTL-Log distribution: 


Vq 


“"V 



2( 

Q—1) 

\/v 

\ 

a+V 

0 ) 



2(«-l) 


1 < a < 2 
a = 1 

0 < a < 1 


(27) 


5.2 rGrL-Geometric distribution 


The pdf and the cdf of the rGTL-Geometric (rGTL-Geo) random variable, respectively, are: 

,, m “ yT~^[^ -{a- 1)(1 - y)Y-Ya - 2{a - 1)(1 - y)] 

J [y-i ^-1 y) n /1/1 (2o) 


F{y,a,u,e) = 


{1 — 6{1 — yy[a — (a — 1)(1 — 

l-{l-yr[a-{a-l){l-y)r 


1 - e{{l - yy[a -(«-!)(!- y)^} 

with 0 < a <2,iy > 0 and 0 < 6^ < 1. After some passages we obtain the hazard rate, given by 

9{y,a,iy) 


(29) 


h{y;a,ye) = 


(30) 


[1 - G{y, a, z/)] ■ {1 - e[l - G{y, a, v)]}' 

The quantile for the rGTL-Geo distribution can be obtained from (16), by considering that A“^[(l — 

i)Am = . 


5.3 rGrL-Poisson distribution 


The pdf and the cdf of the rGTL-Poisson (rGTL-Poi) random variable, respectively, are: 


f{y,a,o,9) = 

Iff) 

^ [« - 2(a - 1)(1 - y)] (1 - yy-^ [a-{a- 1)(1 - y)r^ ■ 
exp{0(l - yY [a-{a- 1)(1 - y)Y} 


(31) 


F(y;a,z/,0) = 1 - 


exp{0(l - yY[a - (1 - a)(l - y)Y} - 1 


exp(0) - 1 

with 0 < a < 2,0 > 0 and 6^ > 0. After some passages we obtain the hazard rate, given by 

9g{y; a, 


h{y;a,iy,9) = 


g6»[l-G(y;a,i/)] _ 


(32) 


(33) 


From the plots reported in Fig. 1, we can state that the hazard rate of the rGTL-Poisson distribution 
can be monotonically increasing, bathtub and UB-BT. 

In this case A~^ [(1 — q)A{9)] = log[g(l — e®) + e% thus the quantile for the rGTL-Poi distribution 
can be obtained by inserting this expression in (16). 
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5.4 rGrL-Binomial distribution 


The pdf and the cdf of the rGrL-Binomial (rGTL-Bin) random variable, respectively, are: 




m9v{l—yY ^[a — {a — 1){1 — yff ^[a — 2{a 

(0 + 1 )™ - 1 

{e{l-yY[a-{a-l){l-y)r + ir-^ 


m-y)] 


( 34 ) 


F{y;a,u,9,m) = 


{9 + 1)™ - {0(1 - yY[a -{a- 1)(1 - y)Y + 1}™ - 1 


(0 + 1)™ 

with 0<q;<2,z/>0, 0>O and m = 1, 2, .... After some passages we obtain the hazard rate, as: 


(35) 


h{y;a,u,9,m) 


m9g{y,a,iy){9[l-G{y,a,iy)] + l}^ ^ 
{0[1 - G{y,a,u)] + 1}™ - 1 


(36) 


Finally, from (16), it is possible to obtain the expression for the quantile of rGTL-Bin distribution, by 
considering that A“^[(l — q)A{9)] = {(1 + g)(0 + 1)™ — g}^/™ — 1. 


6 Applications 

In this section, we fit some models, belonging to the proposed class, to real data. In particular, in the 
first example, we consider the dataset reported in Genf (2013), regarding two different algorithms, 
SC16 and P3, used to estimate unit capacity factors by the electric utility industry, while, in the second 
examples, we consider the percentage of muslim population and the percentage of atheists, used by 
Silva and Barreto-Souza (2014). 


Example 1. In Gen 9 (2013), the author fits the TL distribution to the capacity data. We compare 
the reported results with those obtained considering the rGTL distribution, the Beta distribution and 
three models from the rGTL-PS class. 

In Table 2 the ML estiamates of the parameters, with the corresponding standard errors, and the val¬ 
ues for Akaike Information Criterion (AIC) are reported for both SC 16 and P3 algorithms. The lower 
values of the AIC obtained for rGTL-Log model, compared to those obtained in correspondence with 
all others models, suggest the superiority of the former in describing these data. Furthermore, Table 
2 gives the results obtained from the Kolmogorov-Smimov (KS) test. Once again, the KS statistic for 
the rGTL-Log distribution is the lowest among all those obtained for the considered distributions. 

Finally, in Fig. 2 are shown, for the two algorithms, the fitted density for the rGTL-Log, the TL, the 
rGTL and the Beta model. As it can be seen, also the plots confirm the superiority of the rGTL-Log 
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rGTL-Log 

rGTL-Geo 

rGTL-Poi 

rGTL 

TL 

Beta 




SC16 



a 

1.3980 (0.687) 

0.8856 (0.621) 

0.6184 (0.511) 

0.5444 (0.431) 


a=0.4869 (0.121) 

v 

0.8665 (0.455) 

0.5578 (0.576) 

1.0414(0.544) 

1.5194 (0.518) 

0.5943 (0.1239) 

b=L1679 (0.358) 

e 

0.9920 (0.014) 

0.9055 (0.119) 

2.1089 (1.311) 




AIC 

-16.7599 

-12.6145 

-8.1251 

-8.0792 

-14.2302 

-15.2149 

KS (p-value) 

0.1071 (0.9544) 

0.148 (0.6952) 

0.2376 (0.1491) 

0.3287 (0.0139) 

0.1690 (0.5272) 

0.1836 (0.4202) 


P3 


d 

1.3275 (0.777) 

0.9098 (0.650) 

0.6455 (0.550) 

0.5573 (0.465) 


a=0.5539 (0.142) 

v 

0.9141 (0.475) 

0.6557 (0.583) 

1.0148 (0.562) 

1.4533 (0.523) 

0.6778 (0.145) 

b=L2198 (0.376) 

e 

0.9821 (0.031) 

0.8611 (0.160) 

1.9458 (1.402) 




AIC 

-10.6097 

-8.2475 

-5.3616 

-6.342 

-8.9965 

-9.5638 

KS (p-value) 

0.1345 (0.8212) 

0.1432 (0.758) 

0.2383 (0.1642) 

0.3395 (0.0099) 

0.1848 (0.4400) 

0.2002 (0.3413) 


Table 2: ML estimates of the parameters, AIC values and Kolmogorov-Smimov test results for the 
first example dataset. 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 


Figure 2: Empirical and fitted density functions for SC 16 (left panel) and P3 (rigth panel) algorithms 
data. 


model in both cases. 


Example 2. In this example, we consider the proportions of muslim population in 152 countries and 
the proportion of atheists of 137 countries. These datasets have been considered by Silva and Barreto- 
Souza (2014) with the aim to select the best model between the Beta and the Kumaraswamy. Along 
these two models, we also consider rGTL-Log, rGTL-Geo and rGTL-Poi models. The ML estimates 
of the parameters, with the corresponding standard errors, and the values for AIC are reported in 
Table 3. In both cases, the rGTL-Log model appears to be the best model, as suggested by the lowest 
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0.0 0.2 0.4 0.6 0.8 1.0 


0.0 0.2 0.4 0.6 0.8 


Figure 3: Empirical and fitted density functions for the proportions of muslim population and atheists. 


value for the AIC. We note that, for the Atheism dataset, also the rGTL-Geo seems to have a better 
performance than the Beta and Kumaraswamy distributions. In Fig. 3 the fitted densities for the 
considered models are shown. 



rGTL-Log 

rGTL-Geo 

rGTL-Poi 

Beta 

KW 



Muslim 



a 

1.3997 (0.341) 

0.6521 (0.178) 

0.4573 (0.120) 

a=0.2976 (0.028) 

a=.2715 (0.033) 

z) 

0.2624 (0.063) 

0.1226 (0.100) 

0.5530 (0.086) 

b= 0.5159 (0.058) 

b=0.5906 (0.057) 

e 

0.9997 (3e-04) 

0.9796 (0.018) 

2.5828 (0.448) 



AIC 

-250.545 

-150.499 

-71.596 

-232.908 

-225.667 



Atheism 



a 

0.8411 (0.744) 

0.9952 (0.608) 

0.5353 (0.283) 

a=0.4368 (0.043) 

a=0.5091 (0.042) 

v 

2.3502(1.242) 

0.9730 (0.758) 

3.0655 (0.782) 

b= 3.6347 (0.538) 

b=3.0914 (0.412) 

e 

0.9900 (0.004) 

0.9746 (0.021) 

3.3155 (0.580) 



AIC 

-449.703 

-438.989 

-381.875 

-407.951 

-417.785 


Table 3: ML estimates of the parameters and AIC values for the second example dataset. 
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7 Conclusion 


In many fields of applied science, the observations take values only in a limited range, it is the case, 
for example, of percentages, proportions and fractions. To model this type of data, the statistical 
literature offers very few alternatives, mainly the Beta distribution and only recently some authors 
have recovered the Topp-Leone distribution, proposed in 1955, and the Kumuraswamy distribution, 
introduced in literature in 1980. Certainly, the lack in the literature of distributions with bounded 
support contrasts with the huge presence of distributions with unbounded support. With the aim to 
reduce this gap, in this paper we have proposed a new class of distribution functions with limited 
support, namely rGTL-PS, obtained by compounding the Power Series distributions and the reflected 
Generalized Topp-Leone distribution. The proposed class includes, as special cases, some new distri¬ 
butions with limited support such as the rGTL-Logarithmic, the rGTL-Geometric, the rGTL-Poisson 
and rGTL-Binomial. 

Like the Beta distribution, the shape of the rGTL-PS density function can be constant, increasing, 
decreasing, unimodal and uniantimodal depending on the values of its parameters. Unlike the Beta 
distribution, the hazard function of rGTL-PS is much more flexible since it can be increasing, bathtub 
and N-shape. Moreover, the main advantage with respect to the Beta distribution is represented by the 
fact that the proposed model presents a distribution function in a closed form and the quantiles can 
be easly obtained. Finally, applications to some real data sets highlight the potential of the proposed 
model. 
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8 Appendix 


The partial derivatives of g{yi] a, u) and G{yi; a, u) with respect to a and u are: 

. , , , . ( yAu — 1) 2yi — 1 

gaiyu a, u) = giyu a, u) <; ^^ + 


'a - {a - l){l - yi)] [a - 2{a - l){l - yi) 


Ga{yh a, v) = -vyi{l - yiY[a - (a - 1)(1 - yi) 


v-l 


gu{yi, a, Y = gigju a,v){- + ln(l - yi) + ln[a - (a - 1)(1 - y))] 


Gu{yi;a,u) = -{(!- yi)[a - {a - 1)(1 - ?/*)]}''In {(1 - yi)[a -{a - 1)(1 - yi)]} 


gaaiyi] a, u) = [gaiyi] a, iy)f - g{yi] a, u) 




a - 2{a - 1)(1 - yi)) 


Gaaiyh a, v) = -vY - YvlY - ViYY - (a - 1)(1 - yY 


u-2 


-I . ViY-Y , 2|/*-i 1 , yi9{yi-,a,Y 

gauiyi] a, u) = g{yi; a, u) <( ^^^ ^^ } + 


'a-{a-l){l-yi)] [a - 2{a - l){l - y)] } [a - {a - l){l - yY 


GaYyi',Ci: Y = GYyi',a, Y {- + ln(l - y) + \n[a - (a - 1)(1 - y)] 


guu{yi]a,u) = gYyi]a,Y ^ ^ + ln(l - y) + ln[a - (a - 1)(1 - yi)] \ - 


u 


GuYyi'^oc^Y = Gu{yi;a,u) ln{(l - yi)[a - {a - 1)(1 - y)]} . 
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_ . T-, /-N ^''{d[^-G(yi\a,u)]}A!{e[l-G{yi\a,v)]}-(A!'{e[l-G{yi\a,v)]}\ , „ / .x 

=- (A>|i-G(,.o,.-|||)- -^a„dB,(,) 

the elements of the observed information matrix J{rj) are given by 


a" {e[l-G{yi-,a,u)]} 
A'{e[l-G(yi-,a,i/)]} ’ 


J rv 


d^i{r];y) 

da^ 


-E 


a, v)g{yi] a, v) - [go^hju 

[g{yi-,a,iy)]‘^ 


0^'^Ri{i) Gaiyi-,a,u 


2 = 1 


n 

+ 0 '^ R 2 {i)Gaa{yi] a, v) 
2=1 


Jo 


d^ij-n^y) 

dadu 


-E 


[^(|/^;a,z/)]2 


n n 

Y, Ri{i)G.(v.-, a, a, A + »Y1 aviyii Cl, 

i=l i=l 


Jad 


d‘^i{y, y) 

dado 


^ ^ X] i?i(i)Go(|/i; a, z/) [1 - G(|/j; a, u)] 

2=1 2=1 


J. 


d^i{y,y) 

(9z/2 


E 

2=1 


'Qaaiyi-, Cl, ^^diyi-i Cl, Z^) [^o(l/i, Cl, Z/) 


G^{yi,a,u) + 0^R2{i)G^^{yi]a, 


ly 


2=1 


2 = 1 


Ji/6» — 


d‘^i{y,y) 

dudO 


^2(*) + ^ 61.(1/*; «, Z/) [1 - G(i/*; a, z/)] 

2=1 2=1 


Jee 


d^i{y, y) 


dO^ 


n A'{e)A{e)-[A{e)? 


[1 - G(i/*; a, z/)]6 

2=1 
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